Interplay between edge states and simple bulk defects in graphene nanoribbons 
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We study the interplay between the edge states and a single impurity in a zigzag graphene nanorib- 
bon. We use tight-binding exact diagonalization techniques, as well as density functional theory 
calculations to obtain the eigenvalue spectrum, the eigenfunctions, as well the dependence of the 
local density of states (LDOS) on energy and position. We note that roughly half of the unperturbed 
eigenstates in the spectrum of the finite-size ribbon hybridize with the impurity state, and the corre- 
sponding eigenvalues are shifted with respect to their unperturbed values. The maximum shift and 
hybridization occur for a state whose energy is inverse proportional to the impurity potential; this 
energy is that of the impurity peak in the DOS spectrum. We find that the interference between 
the impurity and the edge gives rise to peculiar modifications of the LDOS of the nanoribbon, in 
particular to oscillations of the edge LDOS. These effects depend on the size of the system, and 
decay with the distance between the edge and the impurity. 



I. INTRODUCTION 

The impact of chemical impurities on graphene 
electronic, optical and mechanical properties has 
been studied extensively, both theoretically 1 and 
experimentally 23 . The effect of impurities has been 
found to be very important, for example a chemical func- 
tionalization of the dangling bonds on the edges, or the 
chemical bonding of add-atoms to some of the bulk car- 
bon atoms may allow to change a semi-metallic or metal- 
lic graphene sample into an insulating one 4 . Moreover, 
scanning tunneling microscopy (STM) images performed 
on disordered samples show large-scale effects around the 
defect centers^, associated with a perturbation of the con- 
ductance of these samples 6 . 

While the majority of theoretical works investigating 
the effects of disorder on the electronic properties of 
graphene focus on transport!^, there have been quite 
a number of theoretical studies focusing on the effect of 
disorder on the local density of states*"^]. Most of these 
works address either the effect of a single impurity in an 
infinite systerrP, or the perturbations of the edge states 
when the impurities are localized directly on the edges^. 
Here we investigate the interplay between the edge states 
and the impurity states. Thus we study the perturbation 
of the edge states in finite size systems (e.g. graphene 
nanoribbons) due to impurities which are situated in the 
bulk. The techniques that we use are tight-binding (TB) 
calculations, the density functional theory (DFT) , as 
well as the T-matrix approximation. These calculations 
allow us to obtain the dependence of the local density of 
states (LDOS) as a function of energy and position, as 
well as the position dependence of the eigenfunctions of 
the system. 

We find that, in the presence of the impurity, the zero 
energy edge state wavefunctions are modified, for exam- 
ple their periodicity along the direction parallel to the 
edge of the ribbon changes. This is responsible for small 
oscillations of the LDOS along the edge, which decay 



with the distance between the impurity and the edge. 
Also, about half the finite energy eigenvalues are shifted, 
and the corresponding eigenfunctions hybridize with the 
impurity state. The largest hybridization and eigenvalue 
shift occurs for a certain state which we denote impurity 
state; the wavefunction corresponding to this state has 
maximal amplitude close to the impurity position. 

Our analysis indicates that a disordered graphene 
nanoribbon (GNR) is a very complex system, for which 
the effect of the impurities cannot be treated perturba- 
tively, nor independently from the rest of the system, and 
for which there is a strong interdependence between the 
edge states and the impurity states. Our observation may 
have important consequences for the analysis of trans- 
port in GNRs: it has been showrP^that the conduction 
in GNRs occurs via the edge states, and thus bulk impu- 
rities may have a strong effect on the edge transport. 

The rest of this paper is organized as follows: in the 
next section we describe the density-functional theory 
(DFT), the tight-binding (TB) and the T-matrix meth- 
ods. Section III is dedicated to a discussion of the results. 
We present the conclusions in section IV. 



II. MODEL AND COMPUTATIONAL 
TECHNIQUES 

We consider a simplified system of a GNR with perfect 
zigzag (zz) edges and a single localized bulk defect (see 
Fig. [I]). The effects of such defect are largest on the edge 
containing sites belonging to the same sublattice as the 
defect site. The positions of atoms in the zz-GNR are 
described by the units vectors a\ and a2 (see Figure 1) 
with ||ai|| = ||a 2 || = 2.461 A, a x = (2.1313,1.2305) A 
and a 2 = (0.0, 2.461) A. The graphene lattice is formed 
by two (A and B) sublattices. We will henceforth take 
the defect to be located on an A atom. We refer to the 
direction parallel to the edge as longitudinal, and to the 
direction perpendicular to the edge as transversal. When 



2 



Longitudinal direction 



Mirror site Edge A 




FIG. 1. Characteristics of a zz-GNR (exemplified for n = 3 
and m = 8) : unit vectors, unit cell, sublattices, directions 
and mirror site definitions. 



analyzing the behavior of the LDOS along the transver- 
sal direction we will consider all the A atoms lying in a 
zigzag path (the red dashed line in Fig. [I]). We define 
as the impurity 'mirror' site the atom closest to the pro- 
jection of the impurity on the A edge; depending on the 
distance between the impurity and the edge, this site can 
fall directly above the impurity, or be displaced laterally 
by one atom (see Fig. [I]); in the latter case there exist 
actually two equivalent mirror sites. 

The zz-GNR is characterized by two chiral numbers 
(n, m). The first chiral number, n, is related to the 
width of the ribbon: 



TABLE I. Ribbon characteristics (chiral numbers, width, 
length, number of atoms), as well as the method of choice 
for each GNR considered. 
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A. Tight-binding method 

For the tight-binding calculations we have performed 
an exact diagonalization of the one-orbital Hamiltonian: 

£ = £^l*>01 + I>i*>oi (5) 

i,3 i 



, x 2 , . .. ,,2 (3n — 1) wo 
«,(„) = -(3n-l)||a 1 .ei|| = -i (1) 

where we have defined the following width constant (that 
will be used in the next section) : 

w = ||ai -ei|| = 2.1313A (2) 

The second chiral number, m, is the number of repe- 
titions along the longitudinal direction and hence it is 
related to the length of the zz-GNR by 

L = m||a 2 -e 2 || (3) 

The number of atoms in a zz-ribbon is dependent on both 
n and m : 

N = An ■ rn (4) 

We present the corresponding characteristics for the set 
of ribbons that we study in this paper in Table 1. While 
we use a fixed ribbon length, we have checked that mod- 
ifying this length does not qualitatively change our re- 
sults. 



where Uj is the hopping parameter between the p z or- 
bitals \i) and \j). The second term is modeling the pres- 
ence of defects on a finite number of graphene sites, via 
the introduction of the on-site potentials V{. In this paper 
we consider the first-order, nearest-neighbor approxima- 
tion (tij is —1 when i and j are nearest-neighbor sites and 
zero otherwise). Taking into account realistic values for 
higher-order hopping processes^ does not change quali- 
tatively the results presented here, though unrealistically 
large values for these hopping terms may strongly affect 
the formation of edge states, as well the interaction be- 
tween the edges and the impurity. These effects will be 
addressed elsewhere^. We assume semi-periodic bound- 
ary conditions (periodic in the longitudinal direction and 
open in the transversal direction). 

When taking into account the effect of a single defect 
localized on the site <i, the Hamiltonian ([5| can be written 
as: 

H = -J2\i)(j\ + V\d)(d\ (6) 

(id) 

where the summation is performed over the nearest- 
neighbor sites. The eigenfunctions of (|6| can be written 
as a linear combination of individual orbit als: 
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i 

ch = (i\k) (8) 

where \k) is the level index corresponding to the energy 
Ek- 

The Hamiltonian ([6| can be diagonalized, and the cor- 
responding eigenvalues E k and eigenvectors (respectively 
Cki) can be obtained numerically. One defines the total 
density of states: 

p(E) = J2S(E-E k ), (9) 

k 

where 5 (E) is the Dirac delta function. While the DOS 
spectrum defined by Eq. ([9| is discrete, for realistic sys- 
tems a phenomenological rounding of the delta function 
needs to be introduced: 

p(E) = ^2f(E-E k ) (10) 

k 

where we take / to be a Lorentzian with a width of 0.05. 

Since we are interested not only in the average den- 
sity of states, but also in the dependence of the density 
of states with position, we define also a local density of 
states (LDOS): 

p(r u E) = ^2\c kl \ 2 f(E-E k ) (11) 

k 

B. Density functional theory 

All DFT calculations presented in this paper have 
been performed using the QuantumWise code which is 
a DFT implementation using pseudopotentials and nu- 
merical atomic orbitals (NAO). For all our calculations 
we have employed the single-zeta polarized basis set and 
the pseudopotentials provided by the package. In or- 
der to find out the defect structure we have performed 
full relaxation calculations: the GNR was allowed to re- 
lax in all directions, and the dynamic convergence crite- 
rion for the relaxation calculation was a maximal force 
of 0.05eV/A, with the energy convergence condition of 
each self-consistent calculation of 10 _5 eV. The carbon 
dangling bonds on the edges have been terminated with 
H atoms, since the dangling bonds are not well treated 
in the DFT framework. We have checked that the stan- 
dard GNR electronic structure (band structure and total 
density of states) is retrieved even in the presence of hy- 
drogenation. 

We have investigated the effect of a vacancy inside rib- 
bons of different widths (n = 5 — 10), but having the same 
length (see Table 1). For the largest ribbon analyzed, we 
have investigated the effects of a vacancy on the edge 
when the position of the vacancy is modified. The differ- 
ent sizes of the nanoribbons used in the calculations im- 
posed different choices for the /c-sampling of the /c-space: 



for greater widths we have employed more ^-points in the 
Monkhorst-Pack scheme (1x3x3 to 1x5x5). The 
sampling has been chosen by performing accuracy tests 
for the calculation of the formation energy of a vacancy 
inside the graphene nanoribbon. Inside a graphene sheet 
a fully relaxed vacancy can adopt two configurations^ a 
configuration in which all the atoms are in the graphitic 
plan (planar configuration) and a second in which one of 
the atoms is slightly out of the graphitic plan (off-plane 
configuration). The second configuration is more stable 
than the first one. The converged values for the forma- 
tion energy of a monovacancy in the first configuration 
inside different nanoribbons were in excellent agreement 
with values obtained in previous works^H. These values 
were not sensitive to modeling the vacancy as a ghost 
atom. 

QuantumWise packages provide visual tools to access 
directly to a broad range of numerical outputs among 
which the most important for the purpose of this inves- 
tigation are the total DOS, the dependence of the DOS 
with energy at a given site, and the dependence of the 
DOS as a function of position at a given energy. 



C. T-Matrix approximation 

The LDOS perturbations in a finite-size system in the 
presence of disorder obtained numerically via DFT and 
TB can be compared to those corresponding to an infi- 
nite system. For such a system, the effect of an impu- 
rity on the LDOS can be evaluated using the T-matrix 
approximation^. This method consists in treating the 
effects of the defect perturbatively, and using diagram- 
matic Green's function techniques. Since this technique 
has been extensively used in the past 9 , the results ob- 
tained via this technique will not be the central focus of 
this paper, but we will use it rather as a cross-check for 
the DFT and TB calculations, as well as to compare the 
LDOS in the presence of a single impurity for the infinite 
and finite-size systems. 



III. RESULTS 

In this section we present the results obtained through 
the TB (section III A), DFT (section III B) and T-matrix 
(section III C) techniques. 



A. Tight-binding calculations 

Here we focus on the results obtained via an exact diag- 
onalization of the TB Hamiltonian. We first classify the 
system eigenfunctions for a defect-free (df) and a defect- 
containing (dc) GNR. Secondly we present the variation 
of the LDOS at zero energy along an axis parallel and re- 
spectively perpendicular to the edge of the GNR. Finally 
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FIG. 2. TB calculated DOS in a df-GNR and dc-GNR (n = 
15). Note the central peak corresponding to the edge states, 
as well as the small low-energy peak corresponding to the 
impurity state. 



we analyze the two-dimensional profile of the LDOS at 
various energies for df-GNRs and dc-GNRs. 



1. The eigenvalue spectrum and the corresponding 
eigenf unctions 

As it has been previously showri^ 19 * 20 * , the DOS spec- 
trum of zz-GNR exhibits a central peak. If in the tight- 
binding model we only take into account the nearest- 
neighbor hopping processes, this peak is located at zero 
energy, and the spectrum is symmetric with respect to 
E = 0. We will focus mostly on this situation in what fol- 
lows. Different values of the higher-order hopping terms 
will be considered elsewhere^. Typical DOS spectra for 
a df-GNR and dc-GNR described by a tight-binding hop- 
ping Hamiltonian with semiperiodic boundary conditions 
(periodic along the longitudinal direction and open along 
the transversal one) are presented in Fig. [9] note the cen- 
tral peak corresponding to the edge states, as well as a 
smaller low-energy peak corresponding to the impurity 
stated 

An analysis of the eigenvalues of the tight-binding 
Hamiltonian reveals that the central peak in the DOS 
corresponds to several low-energy states (see Fig. [3]). 
This analysis also reveals that, in the presence of a defect, 
the eigenvalue spectrum is modified, such that approxi- 
matively half of the finite energy eigenvalues are shifted. 
This shift, which decreases with the value of the impurity 
potential, is maximum for a certain state which we de- 
note as 'defect' state or 'impurity' state (see Fig.[3|. The 
formation of this state at an energy for which there is 
no corresponding state in the clean ribbon is responsible 
for the DOS impurity peak depicted in Fig. [9] As we will 
show in what follows (see Fig.|4]e), the eigenfunction cor- 
responding to this eigenvalue is mostly localized in the 
vicinity of the impurity, justifying the identification of 
this state as 'impurity state'. 

The eigenstates corresponding to the eigenvalue spec- 
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FIG. 3. The lowest-energy eigenstates for a n — 15 GNR 
spectrum for an unperturbed system, as well as for different 
values of the impurity potential. 

trum described above can be classified according to their 
energy and their position dependence (all energies are 
evaluated in units of the hopping parameter t which we 
take to be equal to 1). Thus, for a n = 15 clean nanorib- 
bon, the lowest energy eigenstates range in energy from 
10 -16 or 10 -15 (zero energy in the limit of the numeri- 
cal precision), up to about 10 -2 . While there are always 
only two eigenstates localized exclusively on the edges, 
the number of 'edge' states (states that have a maxi- 
mum intensity on the edge and a non-zero decay in the 
bulk) varies with the size of the system. With increasing 
the ribbon size, the overall number of eigenvalues, and 
consequently the number of eigenstates composing the 
zero-energy peak increases (for narrow zz-GNR, n < 5, 
there are only two edge eigenstates). 

The lowest-energy wavefunctions are localized on the 
edge (see Fig.[4]a,b) and exhibit oscillations along the lon- 
gitudinal direction, with a period inverse proportional to 
their corresponding energy. When the energy increases, 
besides the longitudinal oscillations, the wavefunctions 
also exhibit a non-zero exponential decay in the bulk, 
whose extent is proportional to their energy (see Fig. [4] 
b,c). This is consist ent wit h the traditional edge-state 
picture for graphen d 14 * 19 * 20 *. The highest-energy edge 
states have the slowest bulk decay and very little lon- 
gitudinal variation; we denote these states as bulk-edge 
states. Above a certain energy value (in the considered 
example of the order of 10 _1 ) the character of the wave- 
functions becomes fully bulk, and the direction of their 
modulation changes from being parallel to the edge to 
being transversal (see Fig. [4]d). 

In the presence of the defect, the lowest-energy edge 
states exhibit a modified periodicity with respect to the 
unperturbed case (Fig. [4] f), but no actual hybridiza- 
tion occurs between the defect and the edge states (the 
wavefunctions do not show an increased intensity close 
to the defect position). However, some of the higher- 
energy bulk-edge states and bulk states hybridize with 
the impurity state (Fig. [4] g,h). More precisely, as also 
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FIG. 4. Wavefimction typology as resulting from the TB analysis: (a) Localized zero-energy edge states in a n = 15 df-GNR; 
(b) Low-energy edge states in df-GNR - note the non-zero but small bulk extension; (c) Bulk-edge states in df-GNR - the bulk 
extension is significative; (d) Bulk states in df-GNR; (e) Impurity state for a dc-GNR (V = 10); (f) Perturbed edge states in a 
dc-GNR, note that such states do not hybridize with the impurity, but their edge periodicity is modified; (g) Hybrid bulk-edge 
impurity state for a dc-GNR (V = 10); (h) Hybrid bulk-impurity state for a dc-GNR (V = 2); 



indicated in Fig. |3J we find that the two-fold degeneracy 
of the eigenvalue spectrum is broken by the impurity, 
such that, for each pair of eigenvalues, one eigenvalue 
and its corresponding eigenstates are unchanged, while 
the other eigenvalue is shifted, and the corresponding 



eigenstate exhibits a hybridization between a GNR wave- 
function and the defect state. The defect state, for which 
the wave function amplitude is maximal in the vicinity 
of the impurity is depicted in Fig. |4]e). 



For a high value of the impurity potential, or in the 
presence of a vacancy (which is equivalent to an infi- 
nite repulsive on-site potential) the eigenvalue shift is 
significant only for a small number of states with ener- 
gies close to that of the defect state (see Fig. |3| which 
is approaching zero in this case. For an infinite system 
this translates into a sharp impurity peak in the DOS 
close to zero energ}^. For a smaller value of V, the de- 
fect state occurs at higher energies, and, as one can see 
from Fig. |3j a larger number of eigenvalues are shifted by 
smaller amounts. This corresponds to the formation of 
a broader and less intense peak in the DOS of a infinite 
system in the presence of small V impurities^. 
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We have also observed that the eigenvalue shift and 
the energy of the defect state depend on the distance 
between the edge and the impurity. In Fig. [5] we depict 
the corresponding eigenvalues when the distance between 
the edge and the defect is reduced. 

Note that when the impurity is getting very close to the 
edge, the eigenstate that exhibits the largest eigenvalue 
shift does not correspond to the first negative-energy bulk 
state, but to the first positive-energy bulk state, as de- 
picted in Fig. [5] 



FIG. 5. The lowest-energy eigenvalues for a n = 15 GNR 
and for different impurity positions (here dd- c is the distance 
between the impurity and the ribbons center). 



2. Spatial variations of the zero- energy LDOS along special 
directions 

The characteristcs of the wavefunctions described in 
the previous section have a direct influence on the be- 
havior of the LDOS as a function of position at different 
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energies. In this section we analyze the dependence of 
the zero-energy density of states as a function of posi- 
tion. The effects of the impurity are most pronounced 
on the edge made of atoms of the same sublattice as the 
impurity, while the effects of the impurity on the oppo- 
site edge are much smaller. We focus on the dependence 
of the DOS along one of the two special directions (lon- 
gitudinal and transversal) described in Fig. [I] 

a. The variation of the zero-energy LDOS along the 
edge We should first note that in order to accurately 
study the perturbations induced by the impurity on the 
edge states, one needs to focus on ribbons that are suf- 
ficiently wide (n > 5). From a TB analysis of such dc- 
GNRs we find that the presence of a defect perturbs the 
zero-energy edge density of states such that it oscillates, 
exhibiting minima and maxima. The magnitude of these 
oscillations depends on two factors: the width of the rib- 
bon (as well as the distance from the defect site to the 
edge), and the value of the impurity potential. The form 
of the perturbation depends on the distance from the de- 
fect site to the edge: if this distance is an even multiple 
of wo the edge DOS exhibits a maximum on the defect 
mirror site, while if it is an odd multiple of w^, the edge 
DOS exhibits a minimum (Fig. [6]). This is probably the 
effect of the positioning of the central impurity exactly 
on top of the mirror site for odd chiral number n's, and 
to its lateral displacement by one atom for even n's, as 
described in Fig. [I] As expected, in all ribbons the am- 
plitude of the edge DOS oscillations is proportional to 
the value of the onsite potential. 

The amplitude of the edge state oscillations depends 
also on the distance between the defect and the edge. To 
quantify this effect we plot the dependence of the LDOS 
along the edge for various distances between the edge and 
the defect (see Figure. [7]). While for defects relatively far 
from the edge rather uniform oscillations are observed, 
with an amplitude which is decreasing with the distance 
between the edge and the defect, for defects close to the 
edge, the amplitude of the DOS on the sites that are 
closer to the defect is diminishing drastically with respect 
to the intensity on sites further away. 

The average edge DOS is also affected by the distance 
from the defect to the edge. While, when the defect 
is far away from the edge the mean edge DOS is ba- 
sically unperturbed from the defect-free situation (see 
Fig. [8J, when the distance from the defect to the edge 
is reduced, the mean edge DOS is reduced, as the elec- 
trons are pushed away from the edge by the repulsive 
impurity potential. 

b. The variation of the LDOS along the transversal 
direction In what follows we describe the dependence of 
the LDOS as a function of the distance from the edge, for 
different positions of the defect (see Fig. [l0|). The LDOS 
along the transverse line defined in Fig. [1] that passes 
through the impurity is strongly perturbed, so that it 
exhibits a minimum on the impurity site. Note that the 
LDOS on the impurity mirror site depends strongly on 
the impurity position, decreasing with the distance be- 
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FIG. 6. Dependence of the LDOS on position along the edge 
for ribbons having different zz-GNR widths : (a) odd and (b) 
even n. 



tween the impurity and the edge (see Fig. [To| b). The 
odd/even oscillations in Fig. [To] b) are probably the re- 
sult of the difference in the definition of the impurity 
mirror site for different odd/even distances between the 
impurity and the edge (see Fig. [I]). 



3. Two-dimensional behavior of the LDOS 

We now focus on the two-dimensional position depen- 
dence of the LDOS. In Figs. |11|12| we describe the be- 
havior of the LDOS for a df-GNR as well as for two 
impurity positions and different values of the impurity 
potential. For a clean GNR one observes that the zero- 
energy DOS is confined to edge, and is constant along 
the longitudinal direction (see Fig. [TT^ i) ). The intensity 
of the zero-energy LDOS is decaying drastically with the 
distance from the edge, such that it has a significant in- 
tensity only along the first few atomic rows parallel to 
the edge, while vanishing completely inside the df-GNR. 
In the presence of a defect the zero-energy LDOS is non- 
zero both on the edges, as well as in the vicinity of the 
defect (see Fig. Ill] b) for V = 10). The LDOS amplitude 



close to the impurity is however maximal at a non-zero 
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FIG. 7. Edge DOS oscillations for different distances between 
the defect site and the edge dd- e , for a n = 25 ribbon and a 
V = 10. 



' 0.0025- 



0.0020 • 



CO 

o 

Q 

ra 0.0015- 

0.0010 




Perfect (defect-free) ribbon 
Low potential defect 
High potential defect 



b) 



Mean DOS on all edge atoms 
Edge DOS of an un-perturbed ribbon 




Distance to the ribbon center (in units of w Q ) 



Chiral index, n 



FIG. 9. TB (a) and DFT (b) calculated mean edge DOS as 
a function of the distance between the defect and the ribbon 
center for a n = 15 ribbon (a) and n = 10 (b). 



B. Density functional calculations 



FIG. 8. Mean edge DOS as a function of the ribbon width 
in df-GNR and dc-GNR. The mean value of the edge DOS 
in the presence of the defect is modified only for very narrow 
ribbons. Note the strong decay of the mean edge DOS with 
the ribbon width. 



energy of E « 0.08 for V = 10 and E « 0.34 for V = 2, 
corresponding to the formation of the impurity state de- 
scribed in the previous section. Thus, for V = 10 (Fig. [II] 
c) and E = 0.08 the LDOS shows a maximum of intensity 
on the six atoms in the vicinity of the defect. However, 
for V = 2, at energy E « 0.34 corresponding to the for- 
mation energy of the impurity state, the LDOS exhibits 
a superposition between a bulk GNR state and the defect 
state (see Fig. 
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d) 



In Fig. [T2| we plot the LDOS dependence on position, 
for an energy corresponding to the impurity state, when 
the impurity is close to the edge. Note the reduction 
in the intensity of the LDOS on the edge close to the 
impurity site, as well as the edge LDOS longitudinal os- 
cillations. 



We present here similar types of results as in the pre- 
vious section, obtained this time via DFT calculations. 
Such calculations do not have access to the eigenstates of 
the system, but only to the LDOS. The defect considered 
here is a vacancy; in the tight-binding model this would 
correspond to an infinite onsite potential. 

Similarly to the TB calculations, the DFT calculations 
show that a vacancy perturbs the LDOS on the edge. The 
shape of the edge perturbation changes rapidly with the 
distance to the edge, as we describe in Fig. [13] Thus, 
when the distance between the center of the GNR and 
the defect is an even/odd multiple of wo, the induced 
edge DOS fluctuations have a maximum/minimum on 
the defect mirror site. This pattern breaks down when 
the impurity is very close to the edge, in which situation 
the mirror site always exhibits a minimum of DOS inten- 
sity. The irregularity of the observed modifications stems 
from the smaller size of the system, as well as from the 
calculation errors, inherently larger in the DFT calcula- 
tions. 

We note that the average DOS on the edge decreases 
with the distance between the impurity and the edge (see 
Fig{9]b), in a manner which is qualitatively similar to the 
the one predicted by the TB calculations. When analyz- 
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FIG. 10. a) LDOS along an axis perpendicular to the edges 
of the n = 15 zz-GNR and passing trough the defect, for 
three positions of a V = 10.0 defect (dd- c is the distance 
between the defect and the ribbon center), b) The LDOS on 
the impurity mirror site as a function of the distance between 
the impurity and the edge. 



ing the two-dimensional LDOS profile obtained via the 
DFT methods (see Fig. 14), we note that, similar to the 
TB results, the effect of the impurity consists mainly in 
the modification of the intensity on the six sites surround- 
ing the defect 9 . 





FIG. 11. TB calculated LDOS in n = 15 (a) df-GNR ; (b) dc- 
GNR, for E = 0, for V = 10.0 ; (c) dc-GNR, for E = -0.08, 
V = 10.0 ; (d) dc-GNR,for E = -0.38, V = 2.0. 




FIG. 12. The LDOS as a function of position when the impu- 
rity is approaching the edge. Note the reduction in intensity 
on the edge mirror site. 



C. Comparison between the infinite system and the 
finite size systems via the T-matrix approximation 

To check the consistency of the DFT and TB calcu- 
lations in describing the effects of the impurity on the 
LDOS we compare these results with those obtained for 
an infinite system via a T-matrix calculation 9 . The spa- 
tial dependence of the LDOS close to the impurity for an 
infinite grapheme sheet with a single impurity potential 
is presented in Fig. [15] Note that, close to the impurity, 
the behavior of the LDOS is exactly the same as the one 
predicted by the TB and the DFT calculations. Note 
also that the long- wavelength oscillations that have been 
predicted by analytical calculations to decay as 1/r in a 
regular 2d system, and as 1/r 2 in graphene^, decay too 



fast to be observed in either analytical T-matrix calcu- 
lations, as well as in the TB and DFT calculations. We 
have checked that such oscillations are well captured and 
visible in the T-matrix analysis of a bilayer graphene. 



IV. CONCLUSIONS 

We have studied the effect of a single impurity, mod- 
eled as an onsite potential, on the spectrum and LDOS 
of a zigzag GNR as a function of the impurity position 
and size of the impurity potential. We have found that, 
at zero energy, the amplitude of the DOS on the edge of 
the ribbon diminishes when the distance between the im- 
purity and the edge decreases. Secondly we have noted 
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FIG. 13. The zero-energy edge DOS calculated via DFT for 
various defect-to-center distances, for a n = 10 ribbon. 
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FIG. 14. DFT calculated LDOS plots in a zz-GNR (n = 10): 
(a) df-GNR (E = O.OeV); (b) bulk state (E = -1.0); (c) dc- 
GNR (vacancy state visible at E = 0.0); (d) zoom- in close to 
the defect for the state described in c). 
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